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ABSTRACT 


LOWTRAN IIIB is a FORTRAN computer program for prediction 
of atmospheric optical transmittance, developed at the U.S. Air 
Force Geophysics Laboratory (AFGL). LOWTRAN IIIB was received 
in the modified form developed by Naval Weapon Center China 
Lake for use on the UNIVAC 1110 computer [Ref. 1], and has now 
been interfaced to the IBM 3033 computer. 

Due to compiler storage limitation in the IBM computer the 
atmospheric data are read into common storage at the beginning 
of the program. The two dimenstional block data submodule has 
been replaced with a linear data array, and a new subroutine 
(array) written to reformat the data. The basic logic 
structure is unchanged. 

Comparisons of NPS LOWTRAN IIIB computations with direct 
optical extinctions measurement over Monterey Bay and at San 
Nicolas Island indicate that under these circumstances LOWTRAN 
IIIB considerably underestimates the extinctions due _ to 
aerosols wihile the computed molecular absorption is in 
substantial agreement with experiment. 

LOWTRAN IIIB is presently available as a method for 
predicting atmospheric transmittance at low resolution at NPS 


and is suitable for incorporation in simulations and studies of 


electro-optic weapon/sensor systems performance. 
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I. INTRODUCTION 

LOWTRAN IIIB is a FORTRAN computer program, developed at 
the U.S.Air Porce Geophysics Laboratory (AFGL), which was 
received inthe modified form developsé dy Naval Weapen 
Center. It calculates the transmittance cf the earth's 
atmosphere in the spectral region from 0.25 to 28.57 micrcn- 
eters (350 *«c 40,000 cm-1) at 20 cm-1 spectrai resolution on 
a linear wavenumber scale. Six atmospharic models, which 
are trcpical, midlatitude summer and winter, subaztic summer 
and winter, and the U.S. 1962 standard atmosphere, covering 
seasonal and latitudinal variations from sea level tc 100 kn 
are available, in addition tc a capability cf allowing the 
user to input atmospheric parameters of his own chceosing. 
The program includes four aerosol models which are average 
continental, urban, rural, and maritime, and either hazy (§ 
km visibility) cz clear (23 km visibility) haze conditions 
can ke selected in addition te the capability of selecting a 
particular visibility of the user's chcosing. The mncdel 
acccunts for molecular absorption, molecular scattering, and 
aerosol extinction, plus atmospheric refraction and earth 
curvature effects. 

LCWTEAN ITIB is presently an available method fer pred- 
icting atmospheric transmittance at low resolution at NPS 
and may be incorpcrated in simulaticns and studies of 
electro-cftic weaponysensor systems performance. 

Basically, this thesis describes a program to 

a) Develop a cemputer code to calculate atmospheric 
transpittance for particular wavelengths. 

b) Cevelop a ccde to give transmittance and mclécular 
absorgptance with bandwidth appropriat2 to grey-body sctces. 
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marine crtical path, based 
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II. THEORY OF ATNOSPHERIC TRANSMITTANCE 


Beer's law fer linear propagation of monochromatic radi- 
ation in the linear regime can be expressed as fellows; 


A 
= ] = -y4yyl(y eqn 2.1 
= | (») MIL) (eq ) 
or 
-Mi 
I(z) = 1(0) e (eqn 2.2) 
The ratio itz) is defined as the transmittance fT cf ths path 


Ia) 
cf length z, and u is the total extinction coefficient which 


is the sum of the coefficients for total abscrption and 
non-fcrward scattering 

A =A, + As 
Both the scattering and absorption coefficients can be 
divided intc components due *o the molecules of the air and 
the aercscl particles suspended in it; i.¢. 

Aa= kyt Ka 

Ns = Tm + Ta 
where = molecular absorpticn coefficient 
= aerosol absorption coefficient 
= mclecular scattering coafficient 


ras SS 
0 


= aerosol scattering coefficient 

The relative values cf the four coefficients depend strongly 
on the density and molecular composition of the atmesrhere 
and the ccmposition, number density and size distribution of 
the aeroscls. The order of importance in each of the impor- 
tant atmospheric transmission windows is shown in Table I 
{Ref. 2]. We note from this that scattering ky both 
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mclecules and aerescls is of greater importance in the 


visible, and absorption in she infrared (particularly 8 to 
14 um). Since devices are designed to operates threughcust 
these wirdows, it is important to predict the transmixtance 
of the atmcspherte asa function of wavelength and weather 
conditions. This prediction becomes a compiex ozoblem of 
computer modelling. The problem requires a detinition of 
the ccmpcsition, density and pressure of “*he atmospkeric 
gases, together with the frequencies, line strengths and 
line widths of all the spectroscopic transitions of «hz gas 
molecules and the aerosol constituents, and the number, size 


TABLE I 
Cptical Atmospheric Attenuation Coefficient 


Atmospheric Wavelength for Attenuation Coefficients 
windcw EO systems (jin) in order of importance 
Visible Q.4 - 9.7 Ona, Tm, K 
Near Initrared O.7 = Bag Ga oe Tae Kn 
Middie Infrared eG ssa Ke Gy 2. eas Fm 
Far Infrared 8.0 - 12.0 Ke, Ka» Sa 


and ccempcsition distributions of the particles. 

Over the entire wavelength range from visible to 
infrared the absorption by molacules, scattering by aeresclis 
and abscrption by aerosols are the ‘dominant extinction 
machanisms, and should be considered. In the “windcw" 
tegicns c£ good transmittance, the molecular line absorption 
is relatively small, and may in some regions be ignered. 
The remaining extinction in these cegions is dué to aerosol 
scattering (which varies only slowly with wavelength) and 
"continuum absorrftion" by the molecules. 
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A. SCLECOLAR ABSORPTION 


Te ccmpute the mcnochromatic “«ransmittance of *h¢ atro- 
Sphere, we must first obtain accurate data describing tas 
frequencies, intensity and line shape of all absorrtion 
lines affecting the attenuation. These should be developed 
a ccmpilaticon based on certain constraints dictated by 
<mospheric abundances if we need to develop this sonochro- 
matic capability. 

The approach by AFGL is <o calculate the “ransmittance 
at given wavelength for each transizion having finit= 
absorption at that wavelength for each of the mclecules. 
The summaticn over all the molcules gives a monochromatic 
+transmittance which is appropriate directly for laser prorpa- 
gaticn, but must be degraded by inceéegration overs finite 
bandwidth fcr low resolution predictions suitable to non- 
laser sources. 

To calculate the transmittance T = * aes the absorption 
coefficient K, should be known as a function of frequency 
for each line. 

The fcur 2ssential line parameters for each line are the 
resonant frequency, y, (cm-1), the intensity per absorbing 
molecule, S(cm-1/molecule cm-2), the Lorentz line width 
parameter, {, (cm-1/atm), “he energy of the lower state, 
E"(cm-1),and the line half-with at half maximum, <, which is 
propcrticnal to the fressure. The frequency,jy/, is indepen- 
dent of both temperature and pressure. The moleculas 
absorrpticn coefficient is given by 


a $x 
K, W)= ea oral (eqn 2.3) 
S = [qo dy (eqn 2.4) 
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The pressure broadened line width depends in a comspli- 


cated fashion on temperature; for computation this is 
approximated by the assumption of temperature - indepenient 
collisicn diameters, leading to 


Ph o5 
a = 2, .%) oF (eqn 2.5) 
é 
where P = 1 atm, and 
T = 296 degrees Kelvin. 


The intensity, S, is pressure-independent and its +¢emp- 
erature dependence can be calculated from E" and 


S(T) = S(T.) Qy CTS) , Qr (5) , exp [1.439 en (Ty (2qn 2.6) 
$s 


Q,(T) Ap (T) 

whare Q, and Q, are vibrational and rotaticnal partition 
functions respectively. At lower pressure the collision 
broadening diminishes and Doppler broadening becomes impor- 
tan+. In the intermediate pressure cange the Voigt profile 
obtained by convolving the two profiles is used. Once 
Deppler krceadening dcminates the Gaussian profile may be 
used. Lorentz linewidths ar2 typically .001 to .01 cm-1, 
while the Doppler width may be 0.0003 cm-1. 

The AFCRL calculations [Ref. 14] are based on the 
Absorpticn Line Parameters Compilation, which originally 
listed YB, S(T), A(Ps, Ts) and E" for 130,090 lines in the 
range beyond 1am, for each one of the species water, carbor 
Gioxide, czcne, nitrcus oxide, carbon monoxide, methane and 
oxygen. A new version of the compilation including over 
139,000 lines between 0.2 am and 30 am was recently reported 
by AFGL [Ref. 3]. 
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B. CCNTINOUM ABSORPTION 


Continuum absorption occurs as a result of cclii- 
sional interactions retween molecules; <hat is, csl 
between two H,O molecules and those of other gases (pzinci- 
pally H.C: N, collision, since nitrogen comprises approxi- 
mately 80 % of the air) 

In wavelength regions between the absorption bands some | 
attenuaticn cccurs cf a continuous nature which is attri- 
buted to water vapor. The mechanism of water vapor 
continuum extinction lacks a complete theoretical ¢xplana- 
“ion. At present, it is believed that it results from the 
accumulated attenuaticns of the distant wings of H,O absorp- 
tion lines, occurring principally in the far infrazed pazx 
of the spectrum. Other postulates, Such as that the pheano- 
msnon is caused by other absorption mechanisms invoiving H,9O0 
dimers, remain pcssibilities y2t to be proved. 

Hewever, all that we can do at présent is to acccunt for 

moirically, tas¢d on 
what limited experimental measurements we have <tc qo cn, 


the water vafor continuum phencmenon 


wo 


until better line shape theories become availabie. Ee 
should be e¢mpnhasized that further accurate and well cont- 
rolled measurements are urgently requized in ‘rdér ‘tc 
account fer this phencmenon in teal atmospheric situations 
with confidence. 

A commen formulation used (for example in LOWTRAN III3) 
to account for the water vapor continuum attenuation ata 
fixed temperature has been to define the transmittarce T(y 
as fcllcws; 


as =~ Kiy) x Tange 
T(W =e (eqn 2.7) 


where the attenuation coefficient k(y) is given by 
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K(¥y) = Cy [Pas * GCF -Pa] (2qr 2.3) 


where Pio and P, refer to the water vapor partial pressure 
and agskient pressure respectively (atm), and Gdefines the 
gQuantity cf water vaper per unit pathlength (gm cm-2 km-1). 
The quantities C, and Cc, are generally referred to as the 
self and fereign (nitrogen) broadening cosfficients for 
water vapor. 

Values for C. and Cy/C, have been obtained empirically 
from labcratory measurements. In LOWTRAN versions I through 
III, the guantity C,/C, is assumed to remain constant over a 
given wavelength interval. However, a major additicn in 
LOWTRAN IIIB has been to account for the temperature depen- 
dence of Cs and this will be discussed iz the 8-14 am H,O 
continuum region. The H, O continuum radiacion in the 
3.5-4.2 nm region is of much less importance and wili not be 
futher discussed here. 


1. Temperature Dependence 


The water vapor continuum attenuation coefficient 
has been found to have a significant temperature dependence, 
which was not accounted for in «he previous LOWTRAN computer 
codas. Based 9n the laboratory measurements using sampdles of 
water vapcr at elevated temperatures, an approximate e€apir- 
ical expression was obtained by Roberts ¢% al {Ref. 6] for 
the temperature dependence which is given in Eqn 2.9 helcw. 
It was fcund that the attenuation coefficient due tc the 
water vapor continuum increases as the temperature 
decreases. That is, for a fixed amount of water vapor ina 
given patk, one would expect more absorption at lower *emp- 
eratures and less absorption at higher temperatures. This 
is a scmewhat unusual phenomenon. In practice one finds 
less water vapor in the atmosphere under cold conditions, 
therefore, the effect of temperatura on the attenuaticz in 
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relss, threugh 7h 


the 8-14 ain 
total watez content of the 


region plays two competiag 
path and tne attexsuation coef fi- 
cient. 

The empirical fits tc the wavelength and senperacure 
dependence of the water vapor continuum described in koberts 
et al (Ref. 6] have been in LOWIRAN ITIB with the 
appropriate conversicn of units as follows. 


used 


‘ < A A - a - 
The attenuaticn coefficient C, in gm be ata at 


4 
@ 
r] 
—_ 
& 
= 
B 


296 K is given by tke following expression in *h: 


region: 


C, (v, 296) =4. 184557 8exp(-7. 87x10) 


e . . Q ~\ Pay , 4 
where V is the wavenumber in ca (note that =10 /y1 where 


is the wavelength in wm). 
The temperature dependence of the coefficient C, was 


found to vary as: 


Co(WeT) = Cy (V-296) exp 6.08 Ga -1) (eqn 2.10) 


where T is the temperature in degrees Kelvin. 


2. Nitrogen Broadened Coefficient 


CL./C, in the aktove &quation represents the ratio cf 
n/ -$ 


the fereign (nitrogen) broadening coefficient to the self 
broadening coefficient. 

In LOWTRAN IIIB we 
parameter Cy/C, based on 
Supplement LOWTRAN IIIB (Ref. 7]. 


it is assumed that Cy /Ce 


value of 0.002 for the 
presented by 


use a 
che measurements 


(at 296 K) does not 


no supporting msasurements are 


Here, 
vary with temperature (since 
available). 


18 


A 


v8. gees 


herr’ 


eos 


Re Re 


PS Ea 
Qi 


es 
& "I 
Belt 


Thus, further measurements ar2 needed to determine 


more accurately the magnitude of the parameter Cy /Cs and its 
temperature and wavelength dependence. 


3. Transmission Calculations 


The transmittance due to the water vapor centinuunm 
in the 8-14 wm region is calculated for a norizontal path of 
length RANGE(km) at altitude z using the following ¢xpres- 
sion in LOWTRAN IIIB: 


T(”) = exp[-c (7, 256) xW(z) XRANGE| (eqn 2.11) 


where W(z) is the effective H,0 absorber amoun*+ per unit 
path length (in gm-cm atm km") at alcitude z, and C. (% 
296) is the water varer (self broadened) attenuation coeffi- 
cient obtained from laboratory measurements a+ a temperature 
ef 296 K. 

The quantity W(z) is given by: 


- 29C nf - 
W(z) = wz) [Pao 9£ [8-08 (96,1 40.002 P Pao (eqn 2.12) 
where 
w(Z) = gm cm-2/Km of H20 in the path at temperature T, 


Fuso = H20 partial pressure (atm) at altitude 2, 


P. = ambient (tctal) pressure (atm) at altitude z, and 
T (z) = ambient temperature at altitude z (degrees 
Kelvin). 


Note that the temperature dependence of the attenuation 
coefficiert C,(y,T) given in Eqn 2.10, has been incorporated 
into the expression for Win Eqn 2.12. The reason for this 


is so that the temperature variation over a given atmcs- 
fpheric slant path is weighted equally with the water ccnzent 
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along the path. 

It may be worth contrasting Eqn 2.12 with the 
‘SI correspcnding expression which has be2n used in LOWTPAN [ 
through LCWIRAN III, that is: 


rel ak Seer ped 


w(z) = w(z) Pu +0-005 (7, - Fro eda e185 


C. AEROSCL EXTINCTICN COEFFICIENT & MODELS. 


Scattering of radiation occurs from molecules in the 
s. air, from aerosol particles suspended in the air, and from 
water droplets in feg, cain cr hail. The attenuation cf a 
beam depends on the size and number density distributicn and 
- the refractive index of the particles. Atmospheric *ran- 
a smissicn may be stzcngly influenced or even dominated dy 
_ scattering from aercsol particles. The major parameter 
% determining the interaction is the ratio of particle radius 
ss tc the wavelegth of the radiation. 


1. Agrosol extinction coefficien: 


Aeroscl extinction is sum of the absorption and scattsring 
ms Fy aeroscl particles. Although scattering by aeztosol parti- 
. cles whose radius is smaller than about 0.03 times the 
wavelength of the iight may be calculated by Rayleigh 
theory, for the whcle range of particle sradius Mie thecry 
must ke applied. 
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a. Rayleigh Scattering 


The Rayleigh volume total scattering cceffici¢ent 
may be written as 


ge Scr” (be8? ) 


he INA 6-1P (egn 2.14) 


where n 3; refractive index 
the number density of molecules 


depolarization factor due to scme molecular 


anisctrorgy. 
Following Penndorf {Ref. 8] it can be writsen in 
the fcllowing form which depends on pressuze and temperza- 


ture. 


4,0NT =! 
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Om = 4.R0TXI0) ©7734) CP/ot3) Y Km (eqn 2.15) 
The streng dependence on wavelength means that Rayleigh 
scattering is a small effect for wavelengths longez than the 
visible. 

be. Mie scattering 

This is appropriate for the condition cf particles size 
comparable to wavelength; i.e., for large molecules and 
small droplets. This is the most important mechanism for 


atmospheric scattering. McCartney gives the scattering 
coefficient (Ref. 9]. 


O = Mf he YK ten) a (eqn 2. 16} 
where r 3; the particle radius 
N(x) the total number density of particles in the 


size range dr about cr 
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K (o¢0) the area scattering coefficient, the ratio of 


effective area to geometrical area, which 


depends one = 2nrf and refractive index n. 


2. Ae€excsol Model 


The range of conditions in the boundary layer is 
represented by three diffent aérosol models, which are 
tural, urban or maritime. The first two azr2 grcuped 
together as " averaged continental ". 


a. Maritime 


Maritime aerosol composition and size distribu- 
tions are significantly different from crural and urtan 
aeroscl types. The maritime aerosol component is due =o salt 
particies which are caused py she evaporation of saaspray a 
droplets. The concentration cf particles near the surface is 
strongly dependent on wind speed ( above 7 m/sec ) and the 
size distribution is also dependent on relative humidicy. 
This salt-particle number-density decreases rapidly abcve 
about 50C m [Ref. 10]. 


bt. Rural 


Rural aercsol background is partly +he preduct 
ef reacticns between various gases in the atmosohare and 
partly due to dust particles picked up from “he sc¢aztth 
surface. The particle concentration is largely dependent on 
the histcry of the aicmass carrying th2 aerosol particles. 


c. Urban 


Urban aerosols contain certain additives from 
compusticn produ ts and from industry. Shetzle and Fenn 
assumed an addition of 35% soot-like particits of similiar 
distributicn to the rurai aerosol [Ref. 11]. 


22 


2c ee, ee ee ee ee ee ee 
i a ea 


a 
a : TIT... Oy 
aa A. MCDEL ATMOSPHERES 
. The LCWTRAN code provides a choic2 of six atmespheric 
7 models, Tkese include the 1962 U. S. Standard Atmosphere 
9 pilus five supplementary nodels. Surface level conditicns for 
=. the supplementary mcdels are given in Table II (Ref. 1]. 
a There are also ‘two ccnditional haze mod2is - corresponding 
ten to sé€a level visual ranges of 5 and 23 km - provided as 
ws basic input data for LOWTRAN III. Aezosol attsnuaticn for 
a cther visual ranges is calculated using an interpolaticns 
@ 
tt TABLE Il 
i Surface Lewel Condition for Model Atmospheres 
E Model Latitude Fressure Temp. Aiz Water Ozone 
atmospheres & density vaper 
mont (ab) (K} (3/3) (g/n3) (o/ m3) 
Subartic im 

- wintes 60 N,Jan 1013 257 1372 1.2 4. 15-05 
: SUNMET 60 N,July 1910 287 1220 Gast 4. 95-95 
x Midlacitude 
Ba winter &5 aoe 1018 222 1301 ae 5. 6.02-05 
K summer 45 N,July 1013 294 1191 14.0 6.Cz-05 
£. Trepfic2l 15 N 1013 300 1167 19.0 5.62-05 
.* 
in extzapclaticn precedure which utilizes these two models. In 
a 
s addicicn to the mod2l atmospheres =h2 user has the crtica of 
+: inserting his own model atmcsphere, or of buiiding ancthker 
oY med21 by combining varicus parts of the six s*andazd models. 
= Provisicns are made in the LOWIRAN program for inserting 
® radiosonde data. There are lisitS on the accuracy of 
es LOWTRAN transmittance calculations when «he input is radio- 
: sende data alone. Radiosondes provide verticai profiles of 
; the synoptic meteorciogical parameters of temperature, 
& 23 
.; 
® 


pressure, humidity, and wind. Information on aicrometeoro~ 
logical faréméeters (i.e., aerosol size jistributions, iscal 
oxidant concentrations, etc.) is also needed as input. In 
the akserce of microseteorological input re) 
izts medel1 atmospheres for that ite eet this may not 
necessarily be applicable +o a spacific lecation. FO 
further discussicn or optical paremer 


E. CALCULATICN MODELLING 


In the application it is impossible <tc measure 


the transmittance T (Y) averaged over the spectral ktand- 
width LW). accepted by the receiver, as indicated in che 
equation 
TW = Sire ap (eqn 3.1) 

where is the central frequency in the interval, ay. 
Consequently for many applications one is interesced in 
knowing the ¢transmit*anc2 of the atmosphere averaged over a 
relatively wide spectral interval, «hat is, for low resolu- 
tiON.» 

< Thus the t2rm transmittance is somewhat ambigucus unless 


tion, ay » cover which it is averaged. This is particularily 
*rsue in the case of mclecular absorption, since the absorp- 
tion coefficient Ky 2s a tapidly varying functicn of 
frequency. It is because of the capid variation cf Ky wth 
frequency that che averaged tranmittance T joes nox, ia 
general, obey the simple exponent 


~ \ 
Ty) ~ ay axp[—Km (¥) Am] (eqn 3.2) 
6yv 
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u 
= is qualified by scme indication of <he spectral zesolu- 
Co 


ers see [Ref. 12]. 


at a Single frequency. Instead cn? miasures 


gal law. Tha+ is, 
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where Ky represents “he net mcnochromatic mclecuiar 


absorpticn coefficient. On the other hand the mclecular 
scattering coefficiert (f,) and the aerosoi scat ing and 
absorgticn ccefficients(g, and Ka) are slowly varying func- 
tions of frequency, and the average transmittarce ane ys the 
Simple exporential law provided cnly the direct transmitted 
keam is being observed 

There are four basic appreaches to obtaining a low rese- 
luticn transmittance vaiue for a qivén path through the 
atmosphere due to molecular absorption. These are 

(1) direct measurements over the reguired path, 

(2) measurements in the labcratory under simulated cendi- 
tions, 

(3) lirne-by-line (Monochrematic) calculations -kasedocn 
detailed knowledge cf spectroscopic line parameters which 
are then averaged over the required spectral interval,and 

(4) calculations based on band model techniques (which use 
available laboratory and/or field transmittance measurements 
cr actual line data as a basis). 

From the point of view of computations, methcd 3 
involves a considerable amoun+ of work and computer time, 
and ccnseguently methcd 4 has been used most frequently. 

Basically, it involves a graphical or computerized 
“curve fitting" technique for accessing the stcred transmit- 
tance data for the appropriate absorber and scatterer 
amounts ccmcuted from the input atmospheric path conditions. 
The predetermined variations of transmittance with frequency 
are stcred for the various atmospheric constituents, for 
standard path conditions. 

For a oqaiven set of meteorological conditions and 
selected path, the appropriate absorber and scatterer 
amounts in the required path are computed for each ccompcenen 
and the results used to correct the transmittances fer these 
compcenents. The transmittances for the separate processes 
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(line aksorption, continuum absorption, molectlas scat- 
tering, aercsol scattering and aerosol absorption) are then 
multiplied together to give the overall absorptior. Thar is 


Ty) (total) = T,,(Y) (line absorption) x 

T,p(4 (continuum absorption) x 

T,,(») (Rayleigh) x 

T,, () (aerosol). 

It shculd@ be noted that the LOWIRAN computer ccde is 


desigred tc caiculate the ‘transmittance for spectrai bands. 


It should not be used to calculate transmission ‘fcr laser 
tenes. Suit2bi¢e techniques for computing atmospheric ‘tran- 
smission for laser lines or extremely narrew spectrai hands 


are described in [Ref. 13], [Ref. 14], anc (Ref. 15]. 
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IV. EXAMPLE OF PROGRAM USE 
ie A. PROBLEM 
os As a check cn the performance of the cod¢ a sample calcula- 
ea tion was carried out using the same input data described by 
Ee Shlanta and Cornette (Ref. 1], for which they previde 
ns computed output. 
gi The example is the calculation of the transmittance from 


es 2350 to 2450 cm-1 in steps of 5 ca-1 for a slant path fron 
2.5km to 8.5 km ata zenith angle of 65 degree, for a 
subeartic winter model atmosphere, and a 23 km visual tangs. 
@ The fpregram is initiated with data input on 4 cards 
(lines) defining the conditions of the computation, and user 
oer directicns are given et the beginning of the program listing 
= (Appendix A). 

EB CARD 1 #45 ## 144 

ey 5 ; Sub-artic winter atmosphere 

13; An average cortinental aerosol model 

13; For 23 km visible range 

a 0; Por normal operation, exc. 

: CARD 2 ¥&2#&0KK% HHeEEED 500 HEHE 500****65.000 

e 2; Vertical & slant path between two altitudes 
ie 0; Fer normél operation of the program which 

& selects tke shorter path when applicable 
2.500 3; Okserver altitude (kn) 

by 8.500 ; Scurce altitude (kn) 

as 65.000 ; Zenith angle at H1 (deg.). 

A CARD 3 *#2350.000##2450 .000#####5.000 

= 2350.000 ; Initial frequency (cm-1) 

is 2456.000 ; Final frequency cu 3 

" 5.000 ; Frequency intervals at which 

2° transmittance is printed. 
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CARD 4 *¥*0 
0; Te end data. 


B. OUTPOT FROM LOWTRAN IIIB. 


The cutpfut fer this problem is given in Avpendix C. The 
parameters defining the atmospheric path, model atmcspheczes 
and frequency tange are firs* printed out. Follewing tabu- 
lations aive the aksorber amcunzs for horizontal and 
vertical path. At the heading HORIZONTAL PROFILES there are 
13 columns. The first column gives a running inteaer asso- 
ciated with each level. The second column gives the level 
altituds in ka. The next 8 columns give the equivalent 
absorber amcunts per km fer tke following absorbing species: 
water vapcr, uniformly mixed gas, czone, nitrogen censinuua, 
water vagper continuum (10 ym), molacular scattering, aezosol 
extincticn and UV czone, respectively. The néxt <thres 
columns give the mean refractive index moduius frcm that 
level to the level abceve, the equivalenz absorber amcunts 
per km for the water vapor continuum (4 um) and for nitric 
acid. 

A heading VERTICAL PROFILES is then printed followed by 
15 columns. The first and second columns give the integer 
associated with the levels traversed by the path and tha 
height of the level. Then follow 8 columns which give the 
integrated equivalent absorber amounts from the initial 
altitude to the level above (in the same order as indicated 
above). The next & columns are labelled PSI, PHI, BETA, and 
THETA (see Appendix D). 

The tctal equivalent absorber amounts for each abscrber 
species are then summarized beiow in their appropriates 


units. 
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The seccnd line in the total equivalent abserber amount 
table gives the water vapor continuum amount (4 mm) and che 
nitric acid amount. 

A transmittance table, containing 12 coluans, now 
follows. The first 3 columns give the frequency (cm-1) 
wavelength  @mn), and total transmittance. The next 7 
columns show the individual transmittances due tc water 
vapor, uniformly mixed gases, ozone, nitrogen (4 um) 
continuum, total water vapor continuum, molecular scat- 
tering, and aéeresol extinction. The last 2 columns give 
absorgpticn due to aeéerosols and the cumulative integrated 
absorpticn. The latter quantity can be used +0 deternine 
the average transmittance over any given spectral interval 
within the spectral range covered by the calcula*icn. 
Finally, the total integrated absorption from V1 to V2 is 
printed cut together with the average transmittance cver the 
band. 


29 


ee . 


V. APPLICATION 


The purpese of this implementation of LOWTRAN ITIB is to 
present a simpie methcd of predicting atmospheric <cransmit- 
tance (at lew resolution) which is applicable over a wide 
spectral interval and for a wide range of atmospheric path. 

In this study the bandwidtas of the computations have 
been chosen to match the transmission bands of the filters 
used by th2 atmospheric optics measurement group at NPS 
(Ref. 16}. These filters have been used with grey - bedy 
sources tc give wavelength resolution in a number of trans- 


TABLE III 
Broadband Syster 


center filter grey-body 
ve- randwidth " source detector 
length at 1/2 ht tenp. 
at fia Nee K 
4 019 20619-21053 2800 Si 
0.63 -010 15 748 -16000 " " 
0.84 -010 11834-11976 " " 
1.03 010 9662-9756 “1 uy) 
1.06 -010 9390-9479 e a. 
1.60 098 6064 -6447 " " 
Zed 2097 4549-4759 " Ge 
3.80 400 2500 -2778 1800 InSb, 77K 
3.835 e110 2571-2646 " " 
10.66 2. B90 827-1083 “ HgCdTe, 77K 
11.02 ° 10 879-938 " " 


aittance measurements at Monterey and elsewhere. Pigure 5.1 
shows the passbands cr the grey-body source filters superia- 
a 1000 ft, 
2 


posed on curves of atmospheric transmittance fo 
ita water at 


path at sea level cortaining 5.7 mm of precipit 
79°P, The filter bands are shown as Table III. 


bl 
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Figure 5.1 Bandpass Regions for Grey-body Source. 
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A. LOWTHAN CALCULATIONS FOR THE MONTHLY AVERAGE CATA OF 
MCHTEBEY BAY 


Calculations intended for use in the propagaticn studias 
in che Mecnterey 3ay Area make use of the knowr range of 
conditions and three different relative humidities, 70, 80, 
and 90 %, taking intc account the monthly and yearly aver- 
ages cf aizc pressure, air temperature, dewpoint temperature, 
water vapor density, and zelative humidity. These monthly 
and yearly averages are taken from the <thesis by Guner 
(Ref. 4], and ozone density is translated from that source 
to 6.00E-05 gm/m3, consistent with tha nidlatitude model in 
Table II. In this Table IV the average temperature is used 
which is the mean value cf maximum and minimum average +emp- 
ezatures correspcading to day/night conditions. 
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Table ¥V shows the calculated atmospheric transemictance 
for 3 different cases of relativ2 humidities and t2mpeza- 
tures fcr selected wavebands of the broadband system a+ 
level(1.60, 2.15, 3.80, 3.835, 10.66, 11.02 mm), 13.16 km 
range, and with typical May conditions. Figure 5.2 and 


=e 


w 


Figure 5.3 also represent the atmospheric transmittance of 
Table V. 

By ccmparing the numbers for different relative 
humidity, it can be seen that increasing relative humidity 


TABLE IV 
Menterey Bay Area Average Weather Condition 


neey er nee) = ayes neat oer Tega aes hum. 
mbkars eq. eg. n/t : 
JAN. 1020.0 19°93 6. 12 oad 5} 
FEB. 1020.0 10.57 6.63 7.67 76 
MAR. 1017.5 10.84 6.68 7.67 78 
AFR. 1018.0 11.68 8.34 8.53 77 
MAY. 1017.0 13.34 9.45 9.11 a 
JUN. 1015.0 14.45 11.12 10.16 79 
JUL. 1017.0 14.73 11.68 10.55 83 
AUG. 1017.5 15.29 10.23 10.92 84 
SEP. 1013.5 15.84 14.73 10.55 79 
Cer. WWT.5 14.73 10.01 9.40 75 
NCV. 1019.0 12.79 7.23 7.95 72 
DEC. 1018.0 11.12 6.12 7.37 73 
YEAR. 1017.5 13.07 8.89 8.82 77 
TABLE V 


Ccuparison of Transmittance by R.H & Temperature 


wavenumber Rel. rhe et | tamperature( C) 

An cn-1 70% 80% 90 % 0.0 Ucias} 30.0 
1) 1.69 6064-6447 .823 .823 .823 823 - 823 «223 
2) 2-5 4549-4759 .834 .832 .831 - 833 - 833 - 833 
3) 3.80 2500-2778 .803 .797 .791 2794 799 - 803 
4) 3.835 2571-2646 .813 .808 .803 806 -810 - 814 
5) 10.6€ 827-1083 .836 .815 .793 -801 821 - 839 
6) 11.02 879-938 856 .833 .808 818 - 840 - 860 


decreases the atmospheric transmittance. We can alsc see 
that increasing temperature increases the transmittance; 
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Y significant differences of the transmittance occur due *o ‘ 
i +emperature. The numbers from the 8-14 um region are much 
a ; larger chan the other band. This point is consistent with 
os the thecry cf the temperature dependence from chapter II. 
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Be. OFTICAL TRANSMISSION AT SAN NICOLAS ISLAND, CEWCON-78. 


In May 1978 the NPS Optical Physics and Micrometeorolegy 
Groups participated in a Cooperative Experiment, West Ccast, 
Oceancgraphy and Meteorology, 1978 (CEWCOM-78) in ccofera- 
tion with the measurement program OSP-III run by the USN 
Optical Signatures Pregram on San Nicolas Island. Cp*ical 
transmittance measurements frem the R/V ACANIA +o shore and 
wicrcmetecrclogical reasurements on th2 R/T ACANIA were made 
by NES measurement group, and other meteorological data were 
monitcred in the OSP experiment. Weatner conditions limited 
the data ottained. 

The measured transmittance of the cptical path from the 
R/V ACANIA to site A on the northwest «ip of San Nicclas 


TABLE VI 
San Nicolas Data, CEWCON-78 


Date Time Range (mn Broadband (am 
asc 3.6~4.0 oP Tb 0 
5/15 1126 63.0 %(94.99) 
" 1458 2192 45.2 %(96.14) 
" 16 323 2200 57.8 % (95,31) ; 
" 1649 2260 49.3 %(96.63) 


Island is listed in Tabie VI. The <cransmittances are z=he 


cbserved vaiues, unccrrected for the molecular absozpticn 


measured with filter selection from a gzrey-bedy sourcs. 

The micrometeorclogical data taker by “~he NPS 
Mezecrolcgy Team dusing she CEWCOM-78 experimer+ [Ref. 16] 
and the Daily Statistical Summary of Meteorological Data 
(Ref. 17] are used tc find interpolated temperature, relative 
humidity, and visibility. The numbers in parentheses tepzre- 
sent LOWIRAN predictions of transmittance, including beth 
molecules and aeroscl extinction based on these metecrolo- 
gical input data. 
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The ccmparison is given for the on2 day for which surfi- 


TABLE VII 
Meteorological Data 


Time Brcadband Fress. Temp. Dew T. RH. Water vapor 
R (mbar) (deg.)(deg.) (4%) (g/m3) 

1126 3.6-4.0 1013.2 14.52 10.0 81.0 9.8 

1633 ue ui 14.99 10.9 66.7 

1458 9.0-12.0 " 14.80 10.9 74.67 57 

1649 uy) " 14.93 10.0 65.69 Ss) 


cient mateorological data are available (Table VII). This 
compariscn shows the LOWTRAN predicted transmittance for 
these conditions to be much higher than the measurements. 
The overestimation of transmittance by LOWTRAN suggests shat 
the particle mod2zl ircludad in the code may be inadequate, 
or that the model's total particle number density may be 
underestimated by the code. Calibration errors in the 
measurement equipment have been investigated and are consid- 
ered an unlikely source of the discrepancy. 


C. EXPERIMENTAL SEASORENENT FRONS SARINA TO POINT PINOS 


The current modelling effort was undertaken in support 
of an experimental progam of transmittance measureme 
total extinction on ranges over Monterey Bay from Marina to 
Foint Pincs (June 1979) using laser and broadband (thermal) 
sources. The optical extinction was measured a+ 0.4860, 
0.6328, 1.03, 1.906, and 11.05 am. 

The calculational program has therefore been directed to 
computaticn cf molecular absorption and aerosol extinction 
coefficients for sea level propagation und¢r a range of 
weather ccnditions statistically typical of Monte¢rey. This 
has been dcne for a variety of wavelengths appropriate .to 
the measurements, and averaged over the filter bandwidths 
for the Exoadband sources used (se2 Table IIIL.). 
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The aerosoi particle spectrum was measured on beard the 


ACANIA with a Knollenberg counter. This information was 


Qu 


used to predict a scattszing extinction coefficient, tase 

cn cne peint measurement of the particle size spectrun. 
Fer ccmparison of data it is necessary to subtract «he 

sifects cf molecular absorption from «he total optical 


TABLE VIII 
Extincticn Coefficient from Optical Measurement and LOWTRAN 


dace time optical LOWTRAN act aerescl 
(um) ,' bas oe molecular 4) 
a c a-c 

6/6 1734 1.60 §3) . $F $572 5c (0) 
1714 1.06 21.7 2.68 0.165 21e5 32.4 
1718 1.03 19.4 2.69 2.63 16.8 32.3 
1720 0.84 22.1 Pee Ue) 0.196 21.9 32.8 
1721 0.62 22.5 2.83 0.066 22.4 32 .8 

6/7 14C6 11.02 19.17 1.84 1.52 17.7 762 
1404 10.6 1Z.1 2.01 1.49 10.6 7.72 
1423 1.06 35.5 2.45 Q. 16 35.3 48.7 
1434 4.03 @€3.9 & .87 2.56 41.3 49.5 

S85 17202 20:30 ee deze 18.8 7.83 
1558 10.6 10.6 p Fae. 1.26 9.34 AN Mi) 
1615 1.6 aaa Ze 17 0.07 3562 49.2 
1628 1.06 38.3 2.43 Q.15 38.2 a0 35 
1632 T4623) 35.0 & .69 2.4 325.6 51.4 
1634 0.84% 40.1 Zed2 0.18 39.9 51.0 
1636 O<63 S562 2.45 0.067 4% 3.13 330 


extincticn. This is done using LOWTRAN. The value cf rela- 
tive humidity and water vapor content have been cbtained 
from Frof. Schacher and Dr. Pairall who made the aarcsoi 
measurements [Ref. 18]. . 

The results are presented in Table VIII. The axtinction 
in units of £-02 km-1 as measured optically, *he LOWTRAN 
prediction of cetal extinction, the moleculazs contribution 
computed by LOWTRAN, total measured minus LOWTRAN molecular, 

nd the values calculated from aerosol spectra are shown. 

By ‘he compaziscn sean by che Table VIII, optical 
extinccicn computed from the aerosol spectra is much higher 
+han the LOWTRAN totai extinction, and is very often higher 
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*han the measured optical extinction. This may bd¢ due =o 
the fact that aercscl extinction was not measured over thea 
whole path but at only one point. It was measured cn board 
she ACANIA meteorological mast az a height of approxinstely 
10 m abcve mean sea Jevel, at approximactely the centerpoinc 
of the optical path. The optical path in this ex 

ranged frcem 18.3 m above sea level at Marina ‘to al 
level at the ACANIA, du2 to Earth curvature over 
km path. 

Total aerosol extinction from LOWTRAN is much less *han 
the ex*tinction calculated fron the measured aerosol 
spectrum. Subtracting the calculated molecular component o 
the LCWTRAN from the optically measured extinction yielis 
the aeroscl extincticn ("difi."), which is much araater than 
she LCWIRAN aerosol component but somewhat less ‘than chat 
predicted from the aerosol spectra. 1.2., agreement of 
aeroscl extinction is much better with aerosol spec«rum 
measurement than with the LOWTRAN aerosol extinction. ilies 
appears frotadle that the LOWTRAN model consistently under- 
estimates the severity of aerosol extincticen in the marine 
éenvircnmert. 
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ONCLUSION AND RECOMBENDATION 

The ccmparisons of LOWTRAN IIIB predictions with optical 
Measurements and with computations from particle siz:2 
@istribution measurements show considerable differences. 
The differences could conceivably azsise from measurement 
errors in the optical transmittance end the particle 
measurements; the two physical measurements give cx2sults 
Teasonably close to each other but in both cases differ 
markedly from LOWTRAN IIIB. In view of the consideration 
give to froth measurement techniques it seems most likely 
that LOWIRAN IIIB ccnsistently underestimates ‘the aerosol 
attenuation in the infrared wavebands 

It must first be realized that LOWTPAN IIIB contains 
limitaticns. It does not include the effects of 

1.clcuds,fcg or precipitation 

2.cefraction, scintillaticn or distortion of the prorpa- 
gating beam due to ttrbulence 

3.backgreurd interference frem direct or scattered solar 
radiaticn 

4. surface Level reflections and masking 

5.multifle scattering 

In addition the models of refraction and earth curvature 
are simplified and the atmosphere is considered to te hori- 
zontally homegeneous and stable with no inversion. These 
limitaticns mean that at best LOWTRAN IIIB results are of 
"moderate" accuracy. 

Som2 of these limitations have been removed in the 
continuing ¢vclution cf LOWTRAN. LOWTRAN IV (Feb. 1978) in 
addition to using improved empirical ‘*cansmittance functions 
for gases cf small absorber amount, includes ths effects of 
nitric acid as an abscrber and emitter, and incorpcerates a 
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“radiance mode" to compute the cvadiance from ‘the <¢arth 
surface and the atmosphere. 

LCWTEAN V (Feb, 1980) has been upgraded to include new 
aeroscl mcdels dependent on altitude and relative humidity, 
and a modei for fog attenuation. 

It may be expected that implementation of «hese codes 
(preferbly LOWTRAN VY) may give better agreement between 
experimant and compution and allow better validaticn of 
aeroscl extincticn calculation methods. 

LCWTRAN IV and LCWTRAN V are not currently availabise at 
NPS in card cr tape form, and considerable modification will 
be reguired to adapt them to use on the IBM 3033 system. In 
the absence of later forms of LOWTRAN, the IIIB version may 
still be used to compare transmittances under different 
conditions. In addition, little change has been made in the 
absorpticn in later versions, so that the us¢ cf IIIB to 
compute the atmospheric absorption is valid. This allows 
its use +0 compute aerosol extinction from total cptical 
extincticn by suttraction of the calculated abscrp<on. 

It is recommended that an effort should be made to adapt 
LOWTRAN V for the NPS IBM computer system. 
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DEFINITIONS AND SYMBOLS 


Absorption at frequency V; also average 
transmittance . 

Aeroscl number density 

Equivalent absorber amount per km at level J 
Aftgle of arrival at adjacent level 

Input zenith angle (Cesc ee] 

Angle subtended a+ the earth's center as path 
traverses adjacent levels 

Total angle Subtended path at earth's center 
Equivalent absorber amount per km at level J+1 
CCnversion factor from deqrse to radians 
Wavelengtr dependent ccafficient usad in 
rerractive index expression | J 
Wavelength dependent coefficient used ina 
refractive index expression 

Log absorption coefficient for wate= vapor. 
Log absorption coefficient for uniformly mixed 


as 

fog absorption coefficient for ozone 

AbSorptiocn cofficient for nitrogen ‘ 
Absorption cofficient for water vapor continuum 
Extinction coefficient for molecular scattering 
Extincticn coefficient for aerosol models 
Aerosol aksorpticn, coefficient % gat 
Seee epee es coefficient for ozone (UV and visibie 
te n 

water va pcr amount (pr.cm/km) at level I 

Cew point temperature Cc 

Path length frem level I to level I +1 ; 
Wave number increment at which transmittance? is 
calculated i 

Heigh+ increment from ievel IT to I + 1 
Equivalent absorber amcoun+ per km at pes H1 
E ee ae ae ee amount per km for water vapor 
ry eve 

Equivalent absorber amount per km for catkon 
dioxide, €+c. at level Z(TI) 

etree absorber amount per km for ozcna at 
eve 

Equivaien absorber amount per km for nitrogen 

at level 2Z(I 

Equivalent absorber amount per km for water vapor 
continuum at level Z(TI) 
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